Stripe glasses: self generated randomness in a uniformly frustrated system 
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We show that a system with competing interactions on different length scales, as relevant for 
the formation of stripes in doped Mott insulators, undergoes a self-generated glass transition which 
is caused by the frustrated nature of the interactions and not related to the presence of quenched 
disorder. An exponentially large number of metastable configurations is found, leading to a slow, 
landscape dominated long time relaxation and a break up of the system into a disordered inhomo- 
geneous state. 
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Competing interactions on different length scales are 
able to stabilize mesoscale phase separations and the 
creation of spatial inhomogeneities in a wide variety of 
systems. Examples are stripe formation in doped Mott 
insulators, as found in transition metal oxides (TMO) 
domains in magnetic multilayer compounds ||Q, 
or mesoscopic structures formed by assembling polymers 
in solution and amphiphilcs in water-oil mixtures Js]^. 
In many of these cases the tendency towards a perfectly 
ordered array of domains, stripes etc. is undermined 
by frustrating long range interactions 0. Very often, 
these assemblies exhibit a long time dynamics similar to 
the relaxation seen in glasses. In the context of stripes 
it has been argued that the presence of only very few 
quenched impurities might already cause a strictly disor- 
dered glassy state ||. Furthermore, recent molecular dy- 
namics calculations for charge ordering in TMO found an 
anomalous long time relaxation with a power spectrum 
similar to l//-noisc ||. Indeed, there is experimental ev- 
idence for the formation of intrinsic inhomogeneities and 
even a stripe-glass in high temperature superconductors 
and other transition metal oxides |Ic|-|l7f|. In particular 
slow, activated dynamics as observed in NMR experi- 
ments [ p^|Jl5| l exhibits a striking universality, rather in- 
dependent of the details of added impurities etc. It is 
therefore tempting to speculate that glassiness in these 
systems is self generated and does not rely on the pres- 
ence of quenched disorder, which may of course further 
stabilize a glassy state. 

In this paper we show that the competition of inter- 
actions on different length scales in a uniformly frus- 
trated systems exhibits a self generated glass transition 
due to the emergence of an exponentially large number 
of metastable states. This result is obtained using the 
replica approach of Ref . jl8],|l9| and by solving the corre- 
sponding many body problem using the self consistent 
screening approximation pc|j2l| ]. Since only very few 
examples exists for models which exhibit self generated 



glassiness |^,^| , all these approaches are extremely im- 
portant for a better understanding of glassiness in gen- 
eral. Even though our findings apply to broader class of 
problems than stripes in TMO, we will adopt a language 
which is specific to that problem. ]24j] 

A model for a uniformly frustrated system with com- 
petition on different length scales is given by the Hamil- 
tonian jjj: 



H= l -Ja«x {r ^(x) 2 + (V^(x)) 2 + ^(x) 4 } 



Here, y(x) characterizes charge degrees of freedom, with 
</?(x) > in a hole-rich region, </?(x) < in a hole poor 
region, and <p(x) = if the local density equals the aver- 
aged one. If ro < the system tends to phase separate 
since we have to guarantee charge neutrality, (ip) = 0. 
The coupling constant, Q, is a measure for the strength of 
the Coulomb interaction and characterizes the competi- 
tion between short and long range interactions. In case of 
strongly anisotropic, quasi two-dimensional cuprate su- 
perconductors one expects an anisotropy of the gradient 
term in Eqj|, which we neglect for simplicity. Despite 
the absence of a clean derivation of Eqjl] from the many 
electron Schrodinger equation, we note that it describes, 
on a phenomenological level, many of the major com- 
peting effects which yield in microscopic theories a rich 
phase diagram of inhomogeneous spin and charge struc- 
tures p5p . For Q = and r < we expect at low tem- 
peratures long range ordered charge modulations. As 
shown in Ref. pq |, the Coulomb interaction suppresses 
this ordered state for all Q > and finite T. Instead, 
the system undergoes several crossovers. Most interest- 
ingly, at low temperatures, where |r(T)| < a mean 
field analysis of Eq.nl shows that besides a correlation 
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length, £ = 2 (r + 2 V / C7) , an additional length scale, 
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l m = An (2y/Q — r) emerges, which characterizes the 
spatial modulation of the field correlations ]26|j , where 
r = ro + uT (f 2 )- These modulations are particularly 
relevant for low enough T where r(T) < 0, where mean 
field theory gives locally ordered regions with character- 
istic size l m <C £.We will show that a stripe glass emerges 
in this temperature regime. 

An essential prerequisite for the anomalous dynami- 
cal features of glassiness, like aging, memory effects and 
ergodicity breaking, is most certainly the occurrence of 
a large number of metastable states, 7V ms , separated by 
energy barriers which are large compared to the tempera- 
ture. In viscous liquids undergoing vitrification calorime- 
try suggests N ms oc exp {const.V), were, V is the system 
size. This observation is the heart of an ideal glass tran- 
sition scenario based on random first order transitions 
p7[ which was originally motivated by microscopic stabil- 
ity analyses of structural glasses and mean field theories 
for random Potts-glasses. Below a cross-over tempera- 
ture, Ta, a 'viscous', energy-landscape dominated long 
time relaxation sets in due to the occurrence of exponen- 
tially many metastable states, i.e. the configurational 
entropy, S c — kg logA/" ms , becomes extensive. Because 
of the large barriers between these states, the system 
will get stuck for extremely long times in one of the 
metastable states, i.e. it will freeze into a glass, at some 
temperature Tq < Ta which depends for example on the 
cooling rate. Even though this laboratory glass transi- 
tion is purely dynamical, a key ingredient of the ideal 
glass transition scenario is that the dynamical slowing 
arises from proximity to an underlying phase transition 
at Tk < Tq, where the configurational entropy would 
vanish like S C (T) oc T — Tk ■ If such an ideal transition 
exists, even for an infinitely slow cooling rate freezing 
will occur at Tk since all the liquid degrees of freedom 
die out due to this 'entropy crisis' |p8| . 

Detailed theoretical investigation of this scenario have 
concentrated on systems with quenched randomness. A 
major step forward for studying nonrandom systems was 
made in Ref. |jl^,^9|, where a new replica approach was 
developed. Within this approach, the configurational en- 
tropy for a model of a structural glass without quenched 
disorder was calculated and found to be in good agree- 
ment with computer simulations. fil| 

We will use this approach to determine S c for a sys- 
tem governed by Eq.j]. The key idea is to introduce, in 
analogy to the theory of conventional phase transitions, 
an appropriate symmetry breaking field, ip(r), and to 
compute the partition sum 



Dipe 



-H[ V ]/T-Z f rf 3 2; [^(x)- y (x)] 2 



(2) 



where g — > + . The energy / [tjj] = — T log Z[tp] will 
be low if tjj(r) equals to configurations which locally 
minimize TL. Sampling all configurations of the -0-field, 



weighted with exp ^— / [ip] /t) , is therefore equivalent to 
scanning all metastable states such that 



F = lim 777 



f M exp (-/ [il>] /t) 



(3) 



is a weighted average of the free energy in the 
various metastable configurations, where W = 
J Dtp exp (— g/2 J d 3 xtp 2 (x)) is introduced for proper 
normalization. If there are only few local minima, the 
limit g — > + behaves perturbatively and F equals to the 
free energy, F, of the system. However, in case expo- 
nentially many local minima with large barriers between 
them exist, a nontrivial contributions arises from the 
"0-integral even for g — > + and the averaged free energy, 
F, differs from F. This enables us to identify the config- 
urational entropy, S c , via F — F — TS C . @|l9) For an 
illustration of the corresponding free energy landscape, 
see inset of Fig.l. 

An explicit expression for S c can be obtained within a 
replicated theory flq] with 



T 1 

F (m) = — lim — log — 

5^0 m W 



It follows that F = ^nElnil 

orn 



S c 



m— 1 

1 dF(m 



D ipZ m [ip] . 



which gives: 



T dm 



(4) 



(5) 
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Inserting Z [ip] of Eq. || into Eq. || finally leads to: 
Z(m) = lim J D m ip exp I ~ ]T H [(p a ] /T 

\ a— 1 
m - \ 
E / ^VW(x) , (6) 

a.b—1 J 

with F (m) = — — log Z(m). Eq.^jhas a formal similarity 
to the action of the random field Ising model, obtained 
within the conventional replica approach, which allows 
us to use techniques, developed for this model pl[ . In 
the following we use the self consistent screening approx- 
imation (SCSA) of EqJ| and determine the 
Green's function, Gab (q) = (v? a (q)</' b (— q)), m replica 
space. G a b (q) then determines the partition function, 
Z(m), and correspondingly S c . 

The interaction between different replicas is symmet- 
ric with respect to the replica index suggesting the mean 
field ansatz 



Gab (q) = (S (q)-T(q))6ab + T(q), 



(7) 



with equal diagonal elements, G (q), and equal off- 
diagonal elements, J- (q) |50|.The physical interpretation 
of G (r — r') = (ip(r)ip(r')) as thermodynamic (instanta- 
neous) correlation function is straightforward. On the 
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other hand, T {r — r') = limt-joo ((p(r,t)ip(r' , 0)) can be 
interpreted as measuring long time correlations, arising 
from trapping in metastable minima which, in mean field 
theory, have infinite barriers between them. An analo- 
gous structure in replica space follows for the diagonal 
elements, Eg (q), and off-diagonal elements, E^-(q), of 
the self energy, which are given in the SCSA as: 

E A (q)=2 / _^LzU(pM(p + q) (8) 

with A £ {Q,T}. The screening of the interaction is 
characterized by the collective propagators T>g 1 (p) = 

(uT)- 1 + IT 5 (p) and 2V(p) = ^^(p^p) with 
polarization functions IF4 (p) = J ^ a A (q + p) A (q). 
The set of equations is closed by the Dyson equation: 
G~ x (k) = Gq 1 (k) + Eg (k) for the diagonal elements, 
and T (k) = ~ g g |^j^|^| for the off diagonal elements, 

respectively. Q^ 1 (q) = r + q 2 + Qq~ 2 is the inverse 
Hartree propagator. Note, all momentum integrations 
have to be cut-off at |p| = A, which is of the order of 
an inverse lattice constant. Once the Q a \, and T> a b are 
determined the free energy becomes 

F(m)/(2mT) = trlogCy- 1 + trlogP" 1 - tiEQ . (9) 

After performing the trace in replica space for arbitrary 
integer m and analytical continuation to m — > 1, S c fol- 
lows from Eq.B. One finds immeadiately S c — if f(k) 
vanishes. In the following we discuss the numerical solu- 
tion of this set of coupled integral equations. 
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FIG. 1. Configurational entropy density, S c /V as function 
of T/T K for Q = 0.01 and r = -10 (Tk = 1.7586) and 
Q = 0.001 and r = -6 (T K = 1.0422). In both cases 
u = A = 1 is used. Note the strong Q-dependece of S c - The 
inset shows a typical energy landscape for finite S c . 

In Fig. 1, S c /V is shown for two different Q- values as 
function of T. At Ta the long time correlation function 
T (k) emerges leading to S c > and a glassy dynam- 
ics sets in. The corresponding free energy landscape is 



schematically illustrated in the inset. Just as in mean 
field Potts glasses with quenched randomness p7j , S c 
vanishes at a lower temperature, Tk- At Tk the entropy 
of the amorphous stripe solid equals that of the stripe 
liquid. There is no entropic advantage anymore to be in 
a liquid state, leading to an obligatory glass transition 
no matter how slow the cooling rate. The laboratory 
glass temperature, Tq, will lie somewhere between Tk 
and Ta and cannot be determined within our theory. We 
also find that Tk and Ta are only weakly decreasing for 
increasing Q, see inset of Fig. 2. Both temperatures re- 
main finite for Q — ► 0. However, S c (Q — * 0) — > 0, i.e. 
the fragility oc ^ of the glass vanishes. In other words, 
the larger the modulation length, the smaller is the num- 
ber of metastable states. Due to the 1/r interaction, the 
limit Q — > does not smoothly connect to the behavior 
at Q = 0. If one includes a finite screening length l s 
in Eqjl], we expect that the glassy state disappears for 
Q « / s - 4 - 
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FIG. 2. Momentum dependence of the instantaneous (Q) 
and long-time {J-) charge correlation function for the same 
parameters as in Fig 1 at T = Ta- The inset shows the glass 
liquid phase diagram as function of Q for ro = —10. 

In Fig. 2, the instantaneous (Q (k)) and long-time 
(J- (k)) charge correlation functions, are shown at T — 
Ta- Even though no charge ordering occurs, the pro- 
nounced peaks at finite k demonstrate that there is a 
modulated state with strong short range correlations, 
l m < £. Also, since T (k) < G(k), the modulated state 
exhibits an anomalous dynamics, where long time cor- 
relations are only slightly reduced compared to instan- 
taneous correlations. Introducing a Debye- Waller factor, 
W = — log (J-/Q) for k close to the peak maximum, gives 
W = 0.12 (0.13) for Q = 0.01 (0.001). The modulation 
length (inverse peak positions) is 3.5A -1 (6A _1 ) and the 
correlation length is 45A" 1 (80A" 1 ) for Q = 0.01 (0.001). 

Due to the competing interactions in Eq.j^, an entropy 
crisis occurs, causing a transition into a glass. This 
purely thermodynamic characterization of the spectrum 
of metastable states, is only the first important step for 
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understanding glassiness, and the investigation of dy- 
namical features is an even bigger challenge, because 
it requires going beyond mean field theory. An argu- 
ment based on "entropic droplets" explains quite well 
the phenomenology of viscous liquid dynamics [2^ ] and 
can even be made semi-quantitative |3jJ . Here we ap- 
ply these arguments to the present stripe model. The 
entropic droplet argument recognizes an intrinsic insta- 
bility of the homogeneous metastable solutions, as char- 
acterized by Eq. [7| namely, creating a droplet of one 
metastable solution within another costs free energy that 
can at most scale as a surface energy but the exponential 
number of configurations gives an entropic driving force 
for such a droplet that scales with V. A mosaic structure 
hence will form. The activation free energy of turning 
over a single region can be computed where the entropic 
gain is given by TS C . A renormalization group calcula- 
tion, based on Ref. [R2J , leads to a size dependent surface 



tension <r{R) — uq (i?A) with 6 



reflecting the 



fact that the interface between two states is wetted by 
intermediate states. This analysis leads to an character- 
istic energy barrier AE cx (TS' C (T)) _1 which implies a 
relaxation time obeying a Vogel-Fulcher law 



t cx exp 



DT, 



K 



T-T K 



(10) 



An estimate for the surface tension <jq ~ \ r o\l { U C) f° r t ne 
stripe model yields Rq ps At 1 {V<jq/TS c ) 2 for the droplet 
volume and D = 3Va 2 /(AT 2 T K dS c /dT\ TK ). Using our 
numerical results for S c , this leads to D ss 60 — 200, typi- 
cal for moderately fragile and strong glasses, and droplet 
sizes R w (25 - 50)A" 1 « (5 - 10)Z m . Note that this 
estimate is only qualitative since i?o ~ £ and a real sepa- 
ration of scales never occurs. The droplet picture implies 
that the glass state breaks up into domains of different 
metastable states, separated by wetted surfaces, build by 
intermediate states. This physical picture is very simi- 
lar to the conclusions made in Ref. flQ] based on NMR 
experiments. 

In summary, we have shown that an exponentially large 
number of metastable configurations emerges in a system 
with competing interactions on different length scales, 
leading to a glass transition and anomalous long time 
dynamics. This glass state is self generated, implying 
that the barriers characterizing the activated dynamics 
are rather universal and should not depend on details like 
added impurities but only on the generic interactions on 
short and long scales, i.e. the magnetic exchange inter- 
actions and the Coulomb interaction. Furthermore, we 
showed that the magnitude of the frustration controls the 
fragility of the glass transition. Finally, following Ref. 
p7[ , we argued that the configurational entropy causes 
a break up of the stripe glass into a mosaic of domains 
or droplets, build up by the various metastable states, 
allowing us to estimate time scales of motions. This 



causes an intrinsic inhomogeneity of all relevant corre- 
lation functions, modulation and correlation lengths etc. 
in the amorphous glassy state. 
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